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Purpose: Polymorphisms in the interleukin 1 alpha (ILIA) and IL1B gene regions were previously associated with 
keratoconus in a Korean population. In the present study, we investigated whether the ILIA and IL1B polymorphisms 
are associated with keratoconus in a Japanese population. 

Methods: A total of 169 Japanese patients with keratoconus and 390 Japanese healthy controls were recruited. We 
genotyped one ILIA single nucleotide polymorphism (SNP; rs2071376) and two IL1B SNPs (rsll43627 and rsl6944) to 
compare the frequencies of alleles, genotypes, and haplotypes between cases and controls. 

Results: Statistically significant association was observed for rsl 143627 (—31 T>C) in the IL1B promoter region; the T 
allele of rsl 143627 was associated with an increased risk of keratoconus (p=0.014, corrected p value [pc]=0.043, odds 
ratio=1.38). The C allele of rsl6944 (-511 OT) in the IL1B promoter region had a 1.33-fold increased risk of keratoconus, 
although this increase did not reach statistical significance (p=0.033, pc=0.098). The TT genotype of rsl 143627 was 
weakly associated with an increased risk of keratoconus (p=0.033, pc=0.099, odds ratio=1.52). However, no significant 
differences were found in the allele and genotype frequencies between the cases and controls for rs2071376 in ILIA. 
Regarding haplotypic diversity, the haplotype created by the T allele of rsl 143627 and C allele of rsl6944 was associated 
with a 1.72-fold increased risk of keratoconus (p=4.0*10~ 5 , pc=1.6xl(T 4 ). 

Conclusions: Our results replicate associations reported recently in a Korean population. Thus, IL1B may play an 
important role in the development of keratoconus through genetic polymorphisms. 



Keratoconus (OMIM 148300) is a non-inflammatory 
corneal disorder characterized by progressive thinning of 
the corneal tissue, which can lead to severe visual impair- 
ment. Prevalence rates for keratoconus vary widely across 
geographic areas and studies, ranging from 0.0002% to 
2.34% [1]. The exact etiology of keratoconus remains uncer- 
tain, but the disease is currently thought to be triggered by 
various genetic, as well as environmental, factors. Family and 
twin studies provide convincing evidence of the importance 
of genetic factors in the development of keratoconus. The 
prevalence of keratoconus in relatives of keratoconus cases 
is greater than that of the general population [2-4], and mono- 
zygotic twins show a high concordance rate for keratoconus 
[5-7]. 

To date, many genetic loci associated with keratoconus 
have been reported with linkage analysis: lp36. 23-36. 21, 
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2p24, 3pl4-ql3, 4q, 5ql4-q21, 5q32-q33, 8ql3.1-q21.11, 9q, 
12p, 13q32, 14qll.2, 14q24.3, 15q22.33-q24.2, 16q22.3-q23.1, 
20pll.21, and 20ql2 [8,9]. Of these loci, 20pll.21 harbors 
the visual system homeobox 1 (VSX1) gene. Several studies 
have identified mutations in VSX1 as a cause of keratoconus 
development in multiple ethnic populations [9-17], and VSX1 
could be a strong candidate gene for keratoconus. However, 
a similar number of studies did not find VSX1 mutations to 
be specific to the disease [18-26]; therefore, VSX1 mutations 
may be responsible for a small fraction of keratoconus cases, 
suggesting that other genetic factors have more powerful 
effects on the development of keratoconus. In addition to 
VSX1, several other candidate genes have been reported, 
including superoxide dismutase 1 [16,27,28], lipoxygenase 
[29], transforming growth factor, beta 1 [30], secreted protein 
acidic and rich in cysteine [16], human leukocyte antigen 
[31-33], mitochondrial complex I genes [34], and collagen 
type IV, alpha 3 and collagen type IV, alpha four [35]. In 
recent genome-wide association studies, hepatocyte growth 
factor [36] and 2q21.3 including RAB3 GTPase activating 
protein subunit 1 (catalytic) [37] were reported to be associ- 
ated with keratoconus. However, studies in other populations 
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failed to replicate previously reported associations for some 
of these candidate genes [16,17,21,22,38], and for other candi- 
dates, replication studies have not been reported. 

Interleukin (IL)-l is a proinflammatory cytokine that 
induces the production of cytokines and chemokines and 
plays an important role in inflammatory processes. IL-1 is 
involved in various cellular activities, including cell prolifera- 
tion, differentiation, and apoptosis. The IL-1 family consists 
of two proinflammatory cytokines (IL-la and IL-1J3) and the 
IL-1 receptor antagonist (IL-IRa). These proteins are encoded 
by ILIA, IL1B, and IL1RN, respectively, which comprise a 
cluster spanning 360 kb on chromosome 2ql4. Recently, Kim 
et al. [39] investigated the association of polymorphisms in 
ILIA, IL1B, and IL1RN with keratoconus in a Korean popu- 
lation and reported that IL1B promoter polymorphisms 
rsl 143627 (-31 T>C) and rsl6944 (-511 OT) are significantly 
associated with an increased risk of keratoconus. The study 
also suggested that the heterozygous genotype of the ILIA 
intronic polymorphism rs2071376 (+376 C>A) is associated 
with a decreased risk of the disease. However, no replication 
study has been performed in other ethnic populations. The 
aim of the present study was to investigate whether genetic 
polymorphisms in the IL1A-IL1B region are associated with 
the risk of keratoconus in Japanese patients. 

METHODS 

We recruited 169 unrelated Japanese patients with kera- 
toconus and 390 unrelated healthy Japanese controls at 
Yokohama City University, Ryogoku Eye Clinic, Tokyo 
Dental College Ichikawa General Hospital, and Yamaguchi 
University in Japan. The diagnosis of keratoconus was 
based on slit-lamp biomicroscopic findings in one or both 
eyes, including corneal stromal thinning, Fleischer ring, 
Vogt's striae, or Munson's sign by cornea specialists at each 
institute. The diagnosis was also based on characteristic 
patterns on video keratography: inferior corneal steepening, 
inferocentral corneal thinning, or an asymmetric bowtie with 
skewed radial axis. The mean age of the patients was 33.8±9.6 
years; 76.2% of patients were male, and 23.8% were female. 
Controls were healthy volunteers from a geographic region 
similar to that for the keratoconus patients, and the control 
population was comparable in age and sex to the patients 
(mean age, 33.4±9.4 years; 74.6% male and 25.4% female). 
According to the Online Mendelian Inheritance in Man 
database [40], several diseases have been associated with the 
ILIA and/or IL1B polymorphisms: periodontitis, Alzheimer 
disease, osteomyelitis, end stage renal disease, gastric cancer, 
inflammatory bowel disease, Parkinson disease, and diabetic 
nephropathy. These diseases were not included in the patient 
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and control groups. This study was approved by the ethics 
committee of each participating institute and complied with 
the guidelines of the Declaration of Helsinki. All study details 
were explained to all patients and controls before consent was 
obtained for genetic screening. 

Genomic DNA was extracted from peripheral blood 
samples using the QIAamp DNA Blood Mini Kit (Qiagen, 
Hilden, Germany). Procedures were performed under stan- 
dardized conditions to prevent variation in DNA quality. We 
evaluated three single nucleotide polymorphisms (SNPs) 
reportedly associated with keratoconus in a Korean popula- 
tion: rs2071376 in ILIA, and rsll43627 and rsl6944 in ILIB. 
Genotyping of all SNPs was performed using the TaqMan 
5' exonuclease assay with validated TaqMan primer-probe 
sets supplied by Applied Biosystems (Foster City, CA). Poly- 
merase chain reaction (PCR) was performed in a reaction 
mixture with a total volume of 10 ul containing IX TaqMan 
Universal PCR Master Mix (Applied Biosystems), 24 nm of 
each primer-probe set, and 3 ng genomic DNA. The PCR 
conditions were as follows: 95 °C for 10 min, followed by 
40 cycles of denaturation at 92 °C for 15 s and annealing/ 
extension at 60 °C for 1 min. The probe fluorescence signal 
was detected using the StepOnePlus Real-Time PCR System 
(Applied Biosystems). 

Allele and genotype frequencies were estimated with 
direct counting. Haplotype frequencies, Hardy-Weinberg 
equilibrium (HWE), and linkage disequilibrium (LD) were 
assessed using Haploview 4.1 software [41]. Differences in 
allele, genotype, and haplotype frequencies between cases 
and controls were assessed with % 2 . To obtain a measure of 
significance corrected for multiple testing bias, we used the 
Bonferroni method. A corrected p (pc) value of <0.05 was 
considered statistically significant. 

RESULTS 

Hardy-Weinberg equilibrium tests and haplotype block: The 
three SNPs were in HWE among the cases and the controls. 
Figure 1 shows the strength of LD for the three SNPs in all 
559 participants. ILIB promoter SNPs rsl 143627 and rsl6944 
were located in one haplotype block, and the magnitude of 
LD between the SNPs was extremely high (D'=0.83, r 2 =0.67). 
ILIA SNP rs2071376 was not linked with the ILIB SNPs. 

Allele and genotype frequencies: Table 1 provides the 
genomic locations and allele frequencies for the three SNPs. 
In ILIB, the frequency of the T allele of rsl 143627 was signifi- 
cantly increased in the cases compared to the controls (60.7% 
versus 52.7%, p=0.014, pc=0.043, odds ratio [OR]=1.38). The 
C allele frequency of rsl6944 was also increased in the cases 
compared to the controls (60.7% versus 53.7%, p=0.033, 
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Figure 1. Linkage disequilibrium plot of three interleukin 1 alpha 
to interleukin 1 beta single nucleotide polymorphisms in 559 study 
participants. The D' value and r 2 value (in parentheses) corre- 
sponding to each single nucleotide polymorphism (SNP) pair are 
expressed as a percentage and shown within the respective square. 
Red represents a high-pairwise D' value. 



OR=1.33), although this increase did not reach statistical 
significance after the Bonferroni correction was performed 
(pc=0.098). Conversely, the frequency of the A allele of 
rs2071376 in ILIA was slightly higher among the cases 
compared to the controls, but no significant association was 
observed (27.7% versus 24.9%, p=0.337, OR=1.15). 

The genotype frequencies of the three SNPs are shown 
in Table 2. Individuals carrying two copies of the disease- 
risk allele of rsll43627 and rsl6944 had 1.52- and 1.45-fold 
increased risks of keratoconus, respectively, and a marginal 



significant p value of 0.033 was observed for the TT genotype 
of rsll43627 (not significant after correction). Forrs2071376 
in ILIA, the AA genotype had a 1.76 -fold increased risk of 
keratoconus, but it was not significant. In addition, the AC 
genotype of rs2071376 reported in a previous study [39] was 
not associated with the disease. 

Haplotype analysis: A haplotype consisting of rsll43627 and 
rsl6944 was observed (Table 3). The frequency of haplotype 
TC created by the risk alleles was significantly increased in 
cases compared to controls (60.4% versus 47.1%, p=4.0 x 10~ 5 , 
pc=1.6 x 10" 4 , OR=1.72). In contrast, the frequencies of haplo- 
type CC and TT were decreased in the cases compared to the 
controls with OR=0.03 (CC: 0.2% versus 6.7%, p=2.9* 10 6 , 
pc=1.2xl0- 5 ; TT: 0.2% versus 5.7%, p=2.1* 10" 5 , pc=8.4xl0- 5 ). 

DISCUSSION 

The aim of the present study was to assess whether polymor- 
phisms in the IL1A-IL1B region affect the development of 
keratoconus in a Japanese population. Therefore, we geno- 
typed three polymorphisms in the region. Here, we report a 
significant association between IL1B promoter polymorphism 
rsl 143627 and keratoconus in Japanese patients, suggesting 
that the IL1B promoter polymorphisms contribute to the risk 
of keratoconus in Japanese and Korean populations. 

The SNPs rsl 143627 is at the -31 position in the promoter 
region of IL1B. As promoter sequences are potential sources 
of polymorphisms affecting gene expression, rsl 143627 
may play critical roles in IL1B gene expression, which may 
contribute to the risk of keratoconus. This hypothesis is 
supported by previous studies showing that IL1B promoter 
polymorphisms, especially the T allele of rsl 143627 in the 
TATA box, can enhance the expression of IL1B [42-45]. The 
IL1-B protein has been detected in human corneal epithelial, 
stromal fibroblast, and endothelial cells [46], and expression 



Table 1. Allele frequencies of three SNPs in ILIA and IL1B 



Risk Allele 

C . T _ „. t . P0 "".., „ Gene Alleles Risk Frequency, % OR 

SNPs Chr. tion (Build Gene , . „ „ „ , P Pc ,„_ n/ ^,„ 

3 location (1>2) allele Cases Controls (95%CI) 

(n=169) (n=390) 

rs2071376 £ 1^ 2 113,535,395 ILIA Intron C>A A 27.7 24.9 0.337 , 

OA) (0.86-1.54) 

1 38 

rsll43627 (-31 T>C) 2 113,594,387 IL1B Promoter T>C T 60.7 52.7 0.014 0.043 Q7 _° 

rsl6944 ^ 5 2 113,594,867 IL1B Promoter C>T C 60.7 53.7 0.033 0.098 

1, major allele; 2, minor allele; OR, odds ratio; CI, confidence interval. Allele frequency differences were tested using 2x2 contingency 
tables. 
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Table 2. Genotype frequencies of three SNPs in ILIA and ILIB 



SNPs 


Alleles (1>2) 


Geno- 
types 


Frequency, % 

Cases Controls 
(n=169) (n=390) 


P 


Pc 


OR (95%CI) 


rs2071376 (+376 OA) 


OA 


AA 


8.3 


4.9 


0.117 




1.76 (0.86-3.60) 






AC 


38.7 


40.1 


0.763 




0.94 (0.65-1.37) 






CC 


53.0 


55.0 


0.654 




0.92 (0.64-1.32) 


rsll43627 (-31 T>C) 


T>C 


cc 


15.4 


22.2 


0.064 




0.64 (0.39-1.03) 






CT 


47.9 


50.1 


0.633 




0.92(0.64-1.31) 






TT 


36.7 


27.6 


0.033 


0.099 


1.52 (1.03-2.23) 


rsl6944 (-511 C>T) 


T>C 


CC 


36.7 


28.6 


0.058 




1.45 (0.99-2.12) 






CT 


47.9 


50.3 


0.613 




0.91 (0.63-1.31) 






TT 


15.4 


21.1 


0.115 




0.68 (0.42-1.10) 



1, major allele; 2, minor allele; OR, odds ratio; CI, confidence interval. Individual genotype frequency differences were tested using 2x2 
contingency tables for each genotype by combining the remaining genotypes into one category. 



of this protein has reportedly been enhanced in keratoconus 
corneas compared to normal corneas [47]. In keratoconus, 
keratocyte apoptosis has been suggested to contribute to the 
corneal thinning process [48-50], which suggests that the 
enhanced ILIB expression caused by the promoter polymor- 
phism, rsll43627, can induce the overexpression of IL1-B 
protein, resulting in the increased corneal apoptotic activity 
observed in patients with keratoconus. None of the previous 
linkage studies of keratoconus has ever detected 2ql4 locus 
harboring ILIB [8]. Linkage studies have less power than 
case-control association studies to detect small or modest 
genetic effects [51]. ILIB polymorphisms have conferred 
modest risk for keratoconus (e.g., OR=1.3-1.5) in the present 
and previous studies [39], and therefore, the previous linkage 
studies may not have been able to find 2ql4. 

The previous study showed that only the heterozygous 
genotype of ILIA SNP rs2071376 was associated with kera- 
toconus in the Korean population, whereas the allele frequen- 
cies were not significantly different between the patients and 



Table 3. Haplotype frequencies of rs1143627 and rs16944 in IL1B 




Frequency, % 








Haplotypes 


Cases (n=169) 


Controls 
(n=390) 


P 


Pc 


OR (95%CI) 


rsll43627-rsl6944 












TC 


60.4 


47.1 


4.0xl0- 5 


1.6xl0- 4 


1.72 (1.33-2.23) 


CT 


39.1 


40.6 


0.643 




0.94(0.72-1.22) 


CC 


0.2 


6.7 


2.9xl0~ 6 


1.2xl0- 5 


0.03 (0.003-0.31) 


TT 


0.2 


5.7 


2.1xl0- s 


8.4x10-5 


0.03 (0.003-0.37) 



1, major allele; 2, minor allele; OR, odds ratio; CI, confidence interval. Individual haplotype frequency differences were tested using 2x2 
contingency tables for each haplotype by combining the remaining haplotypes into one category. 



the controls [39]. In the present study, there was no association 
between keratoconus and any allele or genotype of rs2071376 
in the Japanese population. Since rs2071376 was not linked 
with the IL1B SNPs, rsl 143627 and rsl6944, in the Japanese 
and Korean populations [39], the possibility that ILIA SNP 
rs2071376 is primarily associated with the development of 
keratoconus is low. 

Keratoconus is genetically heterogeneous, and detecting 
susceptibility genes could provide useful information 
regarding the etiology of this poorly understood disease. In 
the present study, we found 

that the ILIB promoter polymorphism, rsl 143627, is 
associated with keratoconus in the Japanese population. This 
finding is in line with a previous study in a Korean popu- 
lation, suggesting that the ILIB promoter polymorphism is 
an important risk factor for susceptibility to keratoconus. 
Further genetic and functional studies are needed to clarify 
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the contribution of the IL1B promoter region to the develop- 
ment of keratoconus. 
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